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ABSTRACT 

We derive a simple form for the propagator of a massless, minimally cou- 
pled scalar in a locally de Sitter geometry of arbitrary spacetime dimension. 
We then employ it to compute the fully renormalized stress tensor at one 
and two loop orders for a massless, minimally coupled 0^ theory which is 
released in Bunch-Davies vacuum at t = in co-moving coordinates. In this 
system the uncertainty principle elevates the scalar above the minimum of its 
potential, resulting in a phase of super-acceleration. With the non-derivative 
self-interaction the scalar's breaking of de Sitter invariance becomes observ- 
able. It is also worth noting that the weak energy condition is violated on 
cosmological scales. An interesting subsidiary result is that canceling over- 
lapping divergences in the stress tensor requires a conformal counterterm 
which has no effect on purely scalar diagrams. 
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1 Introduction 



One of the most intriguing and potentially important results of observational 
cosmology is the inference, from Type la supernovae at high redshift [jl], |], 
that the universe is entering a phase of "accelerated expansion" . This phrase 
sometimes gives rise to a misconception that is best explained in the context 
of a homogeneous, isotropic and spatially fiat geometry, 

ds^ = —dt^ + a^{t)dx ■ dx . (1) 

The rate of cosmological expansion is the Hubble constant, 

whereas the deceleration parameter is, 

m -1 - H ■ (3) 

Although the general public sometimes takes "accelerated expansion" to 
mean that H > 0, the actual meaning is rather that g < 0. 

In fact no stable theory can exhibit if > on the classical level. This is a 
simple consequence of the two nontrivial Einstein equations in this geometry, 

3H'^ = SttGp, (4) 
-2H-3H^ = SnGp. (5) 

Here p(t) and p{t) are the energy density and pressure, respectively. Adding 
the two equations gives, 

-2H = 87iG{p + p) . (6) 

The quantity on the right-hand side is non-negative as a consequence of the 
Weak Energy Condition ^ which asserts that no physical observer can see a 
negative energy density. Its mathematical transcription is T^^yW^W^ > for 
all timelike vectors VV^. To derive p + p > for a homogeneous and isotropic 
cosmology, note first that the nonzero elements of the stress-energy tensor 
are Tqq = —pgoo and Tij = pQij. Now apply the condition in the limit that 
becomes null. 
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The Weak Energy Condition not only implies H < 0, it also tells us that 
g > — 1. Although stable theories obey the Weak Energy Condition on the 
classical level it has long been known that quantum effects can give rise to 
violations 0], although those studied so far have been on microscopic scales. 
This is of great interest because the present data is actually consistent with 
g < — 1 0. It is only theoretical prejudice in favor of stability and against 
quantum effects on cosmological scales that results in the usual likelihood 
contours being cut off at g = —1. If the proposed Supernova Pencil Beam 
Survey P| or SNAP experiments were to give a definitive determination 
of g < — 1 we should be forced to abandon this prejudice. The purpose of 
this paper is to give a precise formulation of a field theory in which quantum 
cosmological effects do in fact result in p+p < for an arbitrarily long period 
of time. 

The model has been studied before, as has its potential for violating the 
Weak Energy Condition on cosmological scales [Q. It consists of a massless, 
minimally coupled scalar with a quartic self-interaction. Although gravity is 
not dynamical, neither is the background flat. The scalar is a spectator to 
A-driven inflation. The new feature, which is the subject of this paper, is 
that we can now regulate the model in a simple way that preserves general 
coordinate invariance. The acid test of simplicity is that one can go be- 
yond coincident propagators to explicitly evaluate higher loop graphs which 
involve integrations. We shall demonstrate this by computing all one and 
two loop contributions to the expectation value of the stress-energy tensor in 
the presence of a state which is Bunch-Davies vacuum at t = 0. Our result 
conflrms the conjecture of ref. P] that this model shows super-acceleration, 
i.e., p + p < 0. The model should also be of interest as an exercise in quan- 
tum field theory on curved spacetime in which the computations are pushed 
beyond the one-loop approximation that has been so thoroughly studied p. 

This Introduction is the first of eight sections. In Section 2 we present 
the Lagrangian and describe the various diagrams which comprise the expec- 
tation value of the stress-energy tensor at one and two loop order. Section 
3 gives an explicit and relatively simple expression for the D-dimensional 
propagator, the possession of which is what allows us to apply dimensional 
regularization at arbitrary order. In Section 4 we compute the one loop mass 
counterterm and the one loop expectation value of the stress-energy tensor. 
Section 5 evaluates the leading order, local diagrams that were already ob- 
tained in ref. |^. The power of our dimensional regularization technique is 
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displayed in Section 6 by computing the nonlocal diagrams which could not 
be done previously. Although these terms do not contribute at leading order 
(as conjectured previously 0) they are necessary to make the stress-energy 
tensor conserved and they do probe the nonlocal, ultraviolet finite sector of 
the theory. They also contribute some nonlocal ultraviolet divergences which, 
if uncompensated, would compromise the model's renormalizability. Section 
7 evaluates the expectation value of the conformal counterterm which absorbs 
them. (The need for this was also realized before 0.) The fully renormalized 
result is presented in Section 8. 



2 The model and its stress-energy tensor 



Our model has the following Lagrangian, 



(7) 



The counterterms reside in A£, 



AC 



-6m'(j)'y/^ + 6nR - D{D - l)H')<j)'^/^ 



6A 



--5Zd^(f)di,(j)g 



fJ,U 



-9 



4! 



(8) 



Those on the second line are of order A^, just as in fiat space |T0|, and are 
irrelevant for the purposes of this paper. However, the counterterms on the 
top line are needed to remove divergences at one and two loop order in the 
stress-energy tensor. It turns out that 6m'^ is of order A while 6A is actually 
of order A°. The conformal counterterm, 6^ is needed to remove overlapping 
divergences in the expectation value of T^i, at order A^. Note that this term 
makes no contribution at all to pure scalar interactions on account of the 
fact that R = D{D — 1)H^ in /^-dimensional, locally de Sitter background. 

Only the scalar is quantized. The metric is a non-dynamical background 
which we take to be locally de Sitter geometry with cosmological constant 
A = [D — l)i/^. It is most convenient to work in conformal coordinates for 
which the metric takes the form, 



where 



1 



(9) 
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It is useful to keep in mind the relation to co-moving coordinates, 

r]^-H-^e-"* ^ t ^ H-Hn(fl{rj)) . (10) 

We release the state in Bunch-Davies vacuum at t — 0, corresponding to 
conformal time r) — —H~^. Note that the infinite future corresponds to 
7/ ^ 0~, so the possible variation of causally related conformal coordinates 
in either space or time is at most Aa; = Ar; = H~^. 
The stress-energy tensor is, 



6S, 



matter 



X) 



(11) 



(1 + 5Z) 



4! 



-25^ Ud - 1)//W + 9A<l>'y% - (0'); 



(12) 



We want the one and two loop contributions to the expectation value of this 
operator in the state which is Bunch-Davies vacuum at t — 0. If iA{x; x') 
stands for the bare propagator in this state, the lowest order kinetic energy 
contributions are. 



- -g^^{x)g^-{x) 



dp(j){x)d^(j){x) 



- \g,M9'\x^ {a,a;iA(x;x')L,^^+o(A)} . (13) 



The first term is evaluated in Section 4 and the order A correction is computed 
in Section 6. The lowest order contributions from the potential energy are. 



-9ixv[x)K-\iA[x;x)\ +^ 



iA(a;;a;) +0(A^ 



(14) 



This is evaluated in Section 5. The contribution from the conformal coun- 
terterm is, 

/,2/ 



-25e{(^ - l)^V(^) (^1'/'' 



.X, 
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+ [g,,{x)gp^{x) - bp^b'^)^ (d,d, - rv(x)9,) (n \<P\x)\q)} 

-2Si{iD-l)H^g^,{x)iA{x;x) 



+ 



9,.ix)gP-ix) - 6P^5-,] {d,d^ - rV(x)9„) tA{x; x) + 0{\^)} . (15) 



This term is evaluated in Section 7. 

Obtaining a true expectation value, rather than an in-out matrix ele- 
ment, requires use of the Schwinger-Keldysh formalism ill], Il2[. However, 



this formalism reduces to the usual Feynman diagrams for the lowest order 
graphs. We shall not actually need the full Schwinger-Keldysh formalism 
until discussing order A kinetic energy corrections in Section 5. 

Because the state is homogeneous and isotropic, the nonzero tensor com- 
ponents are (f2|Too|fi) = —pgoo and {fl\Tij\Q) = pgij. This is true as well for 
each of the three contributions just described, so we shall finally report an 
energy density and pressure for each. Although neither the kinetic energy nor 
the potential energy contributions is separately conserved, their sum is. The 
conformal contribution is conserved by itself, as is the cosmological countert- 
erm. However, because these quantities harbor divergences, "conservation" 
must be understood with the ultraviolet regulator still on. That is, we must 
keep the spacetime dimension D arbitrary, which makes conservation read, 
p = —{D — l)H{p + p). This is an important check on accuracy. Another 
important check is the cancellation of overlapping divergences that occurs 
when the various contributions are combined. 



3 The D-dimensional propagator 

The behavior of free massless minimally coupled scalar field has been inves- 
tigated extensively |l^, 0, |l^, |l^, Among the curious properties of 



these particles are the absence of normalizable, de Sitter invariant states [1T3 



and the appearance of acausal infrared singularities when the Bunch-Davies 
vacuum is used with infinite spatial surfaces |T3, |T3[. To regulate this infrared 



problem we work on the manifold T^~^ x R, with the spatial coordinates in 
the finite range, —H~^/2 < x^ < H~^/2. Although the actual propagator 
is a mode sum on this manifold, the small possible variation in conformal 
coordinates renders the first term of the Euler-Maclaurin formula — just the 
integral — an excellent approximation, and the finite spatial range serves 
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merely to cut off what would have been a logarithmic infrared divergence on 
D = 



infinite space. In D = 3 + 1 spacetime dimensions the result is |jT9 

< — ■ \n(ii(x: x')^ + ^ 



iA{x; x) 



D=4 



\2ttJ [yix;x') 2 

where the modified de Sitter length function is,[] 

y{x; x') = Q{ri)Q{ri')H^ \\x — x'\\'^ — {\ri — ri'\ — ie) 



(16) 



(17) 



Neglecting the higher order Euler-Maclaurin terms does not prevent (p^) 
from solving the correct differential equation. The higher terms also drop out 
of quite complicated, nonlinear relations such as the Ward Identity for the 
one loop graviton self-energy . We shall therefore regard the technique as 
valid and confine ourselves to finding the appropriate generalization of (|1^) 
to D spacetime dimensions. 

We seek a function of y{x; x') and the two conformal factors which obeys 
the equation, 



v'"'^^''-\v)^^^ix; x') = i5^{x - x') . 



(18) 



One consequence of de Sitter invariance is that the kinetic operator takes a 
particularly simple form when acting on a function of just y{x; x'). To reach 
this form note first that derivatives of y{x; x') give. 



dy_ 

dx^^ 



Hn{ri) y5^^ + 2n{ri')H/\x^ + 2ieHn{ri')sgn{ri - ri')S^ 



(19) 



d'^y 



dxf^dx'^ 



+4ieHn{T]')sgn{7] - + 4ie5(r/ - V)*^ u*^ 



(20) 



^What is termed "the de Sitter length function" in the literature is, 

x') = 1 — y{x] x') . 

The geodesic length from to x''^, £{x; x'), is related to y{x; x') as follows, 

'1 



yix;x) 



-He{x;x'] 
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Now use Q{r]')HAx^ = 1 - n{r]')/n{T]) to conclude, 



dy dy 
dx^^dx^ 



H^Q^{t]) -y^ + Mj - 4:ieHQ{rj')ysgia{r] - t]') 

H'n'iv)[-2y + A + 2iD-2)^ 

—4ieHfl{rj')sgn{ri — i]') — Aie5{rj — r]'] 



Finally, use the chain rule and substitute (|2ll - ^3D 



d^y 



[22) 



(23) 



(24) 



H'Q''{r^){ {Ay - y')f"{y) + D{2 - y)f{y) - 4^e5(r/ - i)f{y) 



2teHnir]')sgniri - i) 2yf"{y) + Df{y) . (25) 



Before considering possible dependence upon the scale factor let us note 
that the delta function on the right hand side of (|TB|) descends from a factor 
of y^~~ in the limit that e — 0. To see this, suppose f{y) has the form, 



fiy) = -^,+o(y 

y2 



(26) 



where ki is a constant. Acting the kinetic operator gives. 



(27) 



Now multiply this by a test function (of x^), integrate / d^x, and take e — 0. 
For e = 0, terms of order y^~~ diverge at x'^ = x'^, but the singularity is 
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integrable. The only delta function comes from the term proportional to, 
, 1 5{rj — 7]') e ttt 5^{x — x') 



Comparison with (|l^) fixes ki, 



(28) 



^^-{^)(^T^m-l). (29) 



2nJ 

Now consider adding to f{y) a symmetric function of ^1{t]) and ^{rj'). 
The only function which can give the same prefactor of Q^{ri) is a constant 
times the same one that appears in (|16|). The D-dimensional propagator 
must therefore take the form, 

tA{x; x') = f(y{x- x')) + h \Yi{n{r^m)) , (30) 

In the e — > limit the function fiy) must obey, 

H^n^{r^){{Ay - y^)ny) + D{2 - y)fiy) 

-Aie6{r] - v')f'iy) - hiD - 1)} = i6^{x - x') . (31) 

We have seen that getting the correct delta function requires a term of 
the form y^~~ plus less singular powers. Series solution generates the higher 
powers. Defining D = 4 — e gives. 



47r2-f V 2) \y^-i \ 2/ n + f (n + l)!r(2 - f ) 4"+i 

This series solves (|3l|) , for A;2 = 0, but it does not reduce to ([l6|) for e = 0. 
The n = term is not even finite in this limit! The resolution to both 
problems is a series of strictly nonnegative integer powers of y, which cancels 
the divergence and the unwanted terms. This series obeys the homogeneous 
equation up to a y^ term which is canceled by the k2{D — 1) term. The final 
result is. 
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+ 1 



E 

n=l 



r(3 



lT{3 + n-e) 1 
nr(2 + n- f) ~ n + §(n + l)!r(2- 



4"+! 



+ 



1 r(3 - e) 

4rri-i: 



ln(fi(r^)fi(V)) • (33) 



The great advantage of this regularization is that it preserves general 
coordinate invariance (once e is taken to zero). One might think that the 
propagator is unwieldy but this is not so in practice. The really cumbersome 
part is the infinite sum on the second line. But these terms all vanish at 
coincidence {y{x; x) = 0) and they vanish for all y{x; x') ai D = A. So one 
need only retain the higher terms when they multiply something else that 
diverges like 1/e. Note also that one need never worry about large y{x;x') 
on account of causality. 

All valid regularizations must reproduce the result of Vilenkin and Ford 
that the coincidence limit of the propagator contains a finite term which 
grows like ln(f2) 



Ht 16, 17, IS 



To check this note that y{x;x) = at 
coincidence. When a variable vanishes like this in dimensional regularization 
one must always assume e to be large enough that the variable is raised to 
only nonnegative powers. We therefore find, 

H \ ^ ( 1 



lim iA{x] x') 




2e 



-r 3 



(34) 

Note that (0) is exact for arbitrary e. 

For every sort of line in the standard Feynman rules the Schwinger- 
Keldysh formalism has a "+" line and a "— " line of the same sort |[Tl|, |1^ . 
Although vertices involve either all + or all — lines, there are ++, H — , 

— h and propagators. All of them are the same function ( p^ ) of the 

appropriate version of the modified de Sitter length function, 

I X — if IP — (|?7 — VI — "^e)^ 



y++[x;x ) 
y+-{x;x') 
y-+{x;x') 
y—{x;x') 



\x 



f I 



\x — X 



{jj — Tj' + leY 
{rj — Tj' — ieY 
(|?7 - ?7'| + ie) 



(35) 
(36) 
(37) 
(38) 



When the points are spacelike related y{x] x') is positive; when they are 



timelike y{x] x') is negative. So zA++(a;; x') and iA. 



(x: X 



are equal (in the 
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limit that e vanishes) for all spacelike separated points and for 77' > rj. That 
is why the ++ and H — contributions cancel whenever x^' strays outside the 
past lightcone of x'^. Note that inside the past hghtcone the ++ and H — 
propagators are conjugate. 



4 The one loop counterterms 

The one loop (A*^) kinetic energy contribution from ([13|) provides a fine illus- 
tration of how (^) is used to compute dimensionally regulated results. Since 
y{x] x') ~ Q{ri)Q{ri') {x — x'Y vanishes at coincidence, only the n = 1 term 
can contribute after the action of two derivatives,^ 



dpd'^iA{x; x') 



257r2 



H_ 



TjA-e) 
2 - 



(39) 



Contracting with the tensor prefactor in (|T3|) gives the one loop kinetic energy 
contribution, 



2%2 



H_ 



'1 



r)r(4 



2 - 



X] 



(40) 



Although ultraviolet finite, this can be nulled by the following cosmological 
counterterm. 



5K 



47r 



IT 

7^ 



)r(4 



+ 0(A). 



(41) 



Fig. 1 depicts the one loop scalar self mass-squared: 



iM, , (x; x' 

1-loop V ■ 



— iA(x; x) + 8m?' 



5^1 



X — X 



(42) 



It is calculated by using the coincidence limit of the propagator given in Eq. 
(|34D . Because of the finite, time-dependent term we cannot make the self 



^Since the time-ordering commutator term should not really be present we avoid the 
delta function term by first evaluating for x^' ^ and then taking coincidence. Another 
way of getting the same result would be to use the — h propagator. 
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+ M 

Figure 1: The scalar self mass-squared at order A. V denotes the 4-point 
vertex and M stands for the mass counterterm vertex. 

mass-squared vanish for all time. A reasonable renormalization condition is 
that it should be zero at t = 0. This can be enforced by, 

*»^ = -^7r(3-l)+0(A=). (43) 

The resulting renormalized self mass-squared is, 

MLJ^; x') = Htd'^ix - x') . (44) 

The temporal growth of has a straightforward physical interpretation. 
It is only at = that the scalar potential has zero curvature; for a general 
field configuration the curvature goes like 0^. But we know from (|3^) that the 
expectation value of 0^ grows with time as the scalar executes a drunkard's 
walk. Hence the mass-squared must exhibit the same time dependence. 



5 The potential energy contributions 

The two loop contribution from the potential energy is depicted in Fig. 2. A 
simple application of the Feynman rules reveals it to be. 



A 



iA{x\ x] 



+ 



5m? 



iA{x; x) 



(45) 



Substituting our results 
counterterm gives. 



for the coincidence limit and (ESl) for the mass 



Aff^ [ / TT Y r^(3 - f ) 1 



\\Hy 



(46) 
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Figure 2: The potential energy contributions to the scalar stress-energy ten- 
sor at order A. V denotes the 4-point vertex and M stands for the mass 
counterterm vertex. 

Since Too = —pgoo, and T^- = pgij, the energy density p and pressure p are, 

The temporal growth of these terms has the same physical interpretation as 
that of the mass-squared. Since the scalar changes slowly we expect that the 
kinetic energy contributions are subdominant. The next section will confirm 
this expectation. 

6 Order A kinetic contributions 

Fig. 3 depicts the order A contributions from the kinetic energy. It is only 
these graphs which require the full Schwinger-Keldysh formalism, very read- 
able derivations of which exist in the literature [jll], |12[ . The rules themselves 
are simple. For every kind of vertex that might appear in a Feynman dia- 
gram the Schwinger-Keldysh formalism has two kinds of vertices: a -|- vertex 
which is the same as that of Feynman, and a — vertex which is the conju- 
gate. Propagation can take place between any two kinds of vertices using the 
appropriate propagator. Each propagator is the same function of the four 
variables y±±{x; x') given at the end of Section 3. 

Using the Schwinger-Keldysh formalism is straightforward. One simply 
draws the analogous Feynman diagram and then sums over ± variations. If 
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D 



M 



Figure 3: Order A kinetic energy contributions to the scalar stress-energy 
tensor. The derivative vertex is D, V denotes the 4-point vertex and M 
stands for the mass counterterm vertex. 

the operator under study is to be time-ordered (as in our case) then the exter- 
nal lines are +. For the diagrams of Fig. 3 this means that the vertex marked 
D is fixed to be +, however one must sum over both + and — contributions 
in the V and M vertices, of course using the appropriate propagators. The 
result is, 



Since both diagrams of Fig. 3 have the same upper loop, they possess a 
common factor in the first curly bracket. The first term within the final 
curly bracket derives from the left hand diagram, while the second term 
comes from the right hand diagram. 

It is convenient to subsume the complicated, e-dependent constants which 
appear in the propagator (PB|), 




iA{x; x') = a{7 {y{x; x')) + f3 \n{n{7])n{r]')) } . 



(49) 
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Comparison with (^) reveals, 



a 



H 



27V 



H 



4r(2-f) ' 

ir(3 + n-e) 
nr(2 + ra- f)4"+i 



and, 
7(2/) = 



'1 - 



A' 



+ E 

n=l 



1 r(3 + n- 



n + f (n + l)!r(2-f) 4"+i 
In this notation, the terms of the last bracket in Eq. (|48|) become. 



-XiA{x'; x' 



ibn? = -i\ap\n{n{ri')) 



Derivatives of the propagators can be written as, 

dy 



d,zA{x;x') = a{^'{y)-^^+(36p'HQ{v)} 



dxP 

These can be further reduced by noting, 
dy 



dxP 

At length one finds. 



X 



p • 



dpiA{x;x')d„iA(^x;x') = a^^4Q\r])n\r]')H^AxpAx^j 
+n'{r^)H%^S,'\yS' + 2Piy + P 



/2 



(50) 



(51) 



(52) 



(53) 



(54) 



(55) 



Expression (|55|) seems complicated due to the infinite sum in the definition 
of 7(1/). However, we need only retain terms which survive as e — 0, 



lim7'(l/) 



lim(7'(y))^ 



y 

1 

,4-e 



1 



+ 



(2_ 
2?/ 



3-e 



+ 



4y2 



(56) 
(57) 
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One can therefore reduce (155|) to 



,4-e 



+ 



2(2 



yi e y2 



,3-e 



,2-e 



(58) 



Note that we have retained the regularization for terms which can produce 
divergences in T^(a;). 

It is useful to break T^^{x) into a sum of six terms of the general form, 



X / d^x'n'^-'iri') In {n{ri')) F„^"(x, x') . 



The functions F^"{x,x') are defined as follows: 

1 



Fr(x,x') 
Fr{x,x') 

Ffix^x') 

Fr{x,x') 



4:H^n{r]'yAxfAx'' 



yl7 yt-- 



2(2 - -)H^n{7]'yAxPAx'' 



H'^n{r]'fAxPAx'' 



y++ y+- 



Fr{x,x') = S^oS'^ 



2Hn{ri') [SP oAx'' + oAx" 
1 1 



.yl7 yt- 



yl- 



.yl7 yl-- 



(59) 

(60) 
(61) 
(62) 
(63) 
(64) 
(65) 



In the remainder of this section we will evaluate explicitly to illus- 

trate the relevant techniques, and then simply give the final answers for the 
remaining five. 

By removing the factors of Q and H the first integral assumes the form, 



5, "5. - -v.uV'^ 
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X J d''x'\n(n(ri'))n\r)') 



AXpAXa- 

Ax 



8-2t 



AxpAXf, 

Axtl^^ 



(66) 



X 



•'II 2 



Here, Ax^_,_ = 
rj' + ie)"^. Now use the differential identities, 



— ilrj — rj'l — ie) , whereas Ax^_ = \\x — x'\\ — {rj — 



1 



1 



-(4 - 26) 

2(2 - e)2 



1 



-(6-26)^^^ 



A^6-26 ' 

-2e(l-e)- 



A^4-26 



to write, 
AxpAxa 



Aa;8-2^ 8e(l -e)(2-e)(3-e) 



dpd^ 



1 



Ax2-2^ 



(67) 
(68) 
(69) 

(70) 



These results pertain for both the ++ and H — terms. Since the range of x^' 
integration does not depend upon x^ ^ the derivatives can be pulled outside. 



iXH'^-^ r2 2-^T(2 



2%6 



X 



J d''x'\n{n{7]'))n\r]') 



1 -e)(3-e) 
e 



1 



A4t^ 



A4-3 



(71) 



It is at this stage that the parallel treatment of the ++ and H — ends. 
The crucial difference between the two terms is that Ax++ contains an ab- 
solute value, I?] — VI) whereas Aa;^ does not. Double derivatives of x^ i_ can 

therefore result in a temporal delta function. Taking e to zero may or may 
not produce a factor of 5^{x — x'), depending upon how many powers of a;+_|_ 
there are. The unique power for general dimension D — A: — e turns out to 
be. 



Axl-^ 



2e(2 - t)i6{A7]) 
\\x — x'W + e2)2~f 



4i- 



7r 



8''{x-x') 



In contrast the H — term gives zero. 



AxX-I 



. 



(72) 



(73) 
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By using ( [73D we can write the H — term in ( [7I| ) in a form that remains 
finite in the unregulated hmit, 





1 


^2 " 


e 




e 



1 



/i- 



2-2e 



Ax 



2-e 



(74) 



The arbitrary mass parameter has been introduced to maintain the correct 
dimensions. Since this term is manifestly ultraviolet finite we may as well 
take e to zero, 



e 



Ax 



2-2e 



In (/x^Ax? 

2Axi_ 



(75) 



The analogous reduction of the ++ term gives, 



e 



In (yU^Axi 



+ ^!!r'^.7 ^^(x-xo, (76) 



er(i-t) 



2Ax^ 



(77) 



Using ([7T|), ( [75| ) and ( [77D we can bring the first integral to the form, 
tXH^-^n-^+'ir]) r(2 - e)T\2 - f ) 



2% 



X / rfVln(^](r/'))fi(^') 1^' 



l-e)(3-6) 

In (^/i^Ax 



- r]i,^d'^ 
In (/i^A, 



2 

++ 



2Ax? 



2Ax? 



era 



r(x-x') . 



(78) 



Now note that the integral has no dependence upon x, so only the temporal 
derivatives matter. Evaluating (5-function integral and taking e to zero in the 
finite terms gives, 

XH^-^Q^+^iv) (1 - f )r(i - 6)r(2 - f )/i- 



2H 



A-e 



e(3 



[ In (^/i^ Ax^_^) In (/i^Ax^ 



Ax? 



Ax^ 



(79) 
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To evaluate the integral in Eq. ([79| ) we pull out another derivative using 
the identities: 



dHn^{Ax^) = 8 
aMn(Aa;2) -- 



In (Ax^ 
Ax2 



Ax2 



Ax2 



(80) 
(81) 



The remaining integrand possesses only logarithmic singularities. If we define 
the coordinate separations, 



Arj = rj — 1]' , r = \\x — x'\\ , 

the logarithms can be expanded as, 

In L2Ax^_^j = In [/i^(Ar/^ - r^)] + me{/\rf - r^) 
In [/i^Ax^_j = In [/i^(Ar7^ - r^)] - m6'(Ar7^ - r^) 

Putting it all together gives the following reduction. 



(82) 

(83) 
(84) 



r In (/i^ Aa;^+) In (/x^/^^ 



Ax? 



Axj 



-il-K'dl I ^ dr]' In {n{r]')) n^{7]') j drr^ (in [jJ^Ar]^ - r^) 

dr]'\n{n{r]'))n^{r]'){A7])^ 
-tSn'^d^ J\ dr]' In {Q{r]'))Q'^{7]') In {2fiAr]) 



2, ^ X 11 
-ln(2/iAr7)- — 



(85) 



The lower limit of temporal integration at rj' = —H^^ (that is, t' = 0) 
derives from the fact that we release the state in free Bunch-Davies vacuum 
at this instant. Ji can now be recast as: 



e(3-e) 



3ln{n{r])) + 



g^2^4(^) ^g J\ In mv')) nivT In {2^iAij) | , (86) 
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where ( is defined as 



TT 



i-o r(i-e)r 1-- 



(87) 



To evaluate the integral in Eq. (p6|), we change variables to ^{r]') = —j^ 
and expand the logarithm, 



1 , 



2 r 



y — 



ln(^)-ln(l]) 



rfl]'^'" In (fi') 



-ln2(n) - -ln(fi) + ln 



2^ 



(n + iy 



(88) 



(89) 



Our final result for Ij^i^if]) is 



^ + lnf^ 



ln((]) 



c 



+ i _ ![! + i V (^-l)(^-2) n-n-l' 
n « « ' ■ ' " 

n=l 



9 6 ■ 6±:^ (n + l)2 



(90) 



The five other terms in Eq. (^) can be evaluated similarly. The results are, 



r4-e 



!^fi-i+T(2-|)r(3-e)r(3-| 



1 



X / rf^x'ln(^](r/'))fi (V) 



In(fi) 



5 7r2 



OO Q— n— 1 




5m '^^^ + 2^^*^ 



4 6 



In(fi) [ . 



(91) 
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ilM 



-T^i2- -)V{Z-e) 



X 



AXpAXcr AXpAXcr 



A4+ 



AxX_ 



2H^'W\l>-m-^,^l^r'-l,lA. (92) 



.!±[L_l]-i+T2(2-i)r(3-6) 



X j d^x' In (^(77')) n\rj') [5p °Ax^ + 5^ °Axp 



2 In (Q) 



3 tiin + iy 



(93) 



3e 



287r6-f6 



r2(2- J)r(3-e) 



X 



J d^x' In [5p °Ax<, + 5^ °Ax 



AxX+ Ax\_ 



2%4 



In (^^) - ^ + ^^^-^ 



(94) 



2%' 



X /d^x'ln(Q(r7'))^l^(V)|^ 



Axt^' Axtl'^ 



+ O^/^^' 



ln(Q) 



00 n— 1 



^ ""6" + ;^,(;^, 



(95) 



Summing the six terms gives the total order A contribution from the 
kinetic energy, 



5 4-'") 
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3 9 



1 ~ (n + 2) ^_^_. 



TT 



oo 



1 



E 



in 



4) 



9 12 e^t'Jn + ip 
We therefore obtain the following order A kinetic energy density, 



PK 



26^ 



and pressure, 
Pk 



267r4 



7 11 



12 " ^ 



9 



2 + 1)^ 



+ In 



2/i 



In (^]) 



ln(r]) 

(^2 - 4) 



7r 
12 



6„^Jn + l)2 



(96) 



(97) 



(98) 



7 The conformal counterterm 



Since C = 1 + 0(e), contains a divergence proportional to gi^,y\\i{yL) / e. 
This is an overlapping divergence, and it must be canceled if the stress- 
energy tensor is to be a well defined operator. It cannot be absorbed into 
a renormalization of the cosmological constant on account of the factor of 
ln(f2). The difficulty is resolved by the conformal counterterm (p!5D,p| 



+ 



-25i{{D-l)H^g^,{x)iA{x-x) 

'g,,{x)g^^-{x) - 5P^5\] {d,d„ - rV(a;)9„) t^{x- x)} 



(99) 



(See the first of the graphs on Fig. 4.) Recall that this comes from a term in 
the action which is contrived to vanish in de Sitter background so it affects 
only the stress tensor and makes no contribution to purely scalar processes. 
The affine connection for de Sitter conformal coordinates is, 



pO 



(100) 



^The potential problem and its resolution were adumbrated in previous work 
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Our simple, dimensional result for the coincidence limit ( p4D of the prop- 
agator implies the following result for the covariant derivative, 



H 



4-e 



"JC • "lie 1 



Substitution and a few trivial manipulations yield 



r(3-e) (?^.+5/5/fi^(r/) . (101) 



4-e 



227r2-e/2 



^)r(3-i)+r(4-e) 



+r(3-e) [g^^u + 6^%,'Q\ 

It is straightforward to extract the energy density and pressure, 
6^H^-' f /3-e 



(102) 



PC 

Pc 



22 ^2-e/2 ^ ^ g 

6^H^-' f /3-e 

~22 7r2-e/2' 



r 3 - 



+ ^(4-e)(ln(^]) 



(103) 



r 3 



+r(4 - e) (In (fi)) - 1 ) + r(3 - e) .(104) 



8 The fully renormalized result 



Summing the potential and kinetic energy densities (|i7|j97|) gives 
XHU /7r\^r2(3- 



PPK 



TT 



+ 



^ + lnf^ 



^-^ln(^^(r])) 



+ i In^ (fi(r/)) + ^ In (r](r^)) - y + ^ + ^^(r/)^' 
_1 - jn + 2^ 



(105) 



Doing the same for the potential and kinetic pressures (|47| , |98|) yields 



Ppk 
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TT- I Vif2 



r2(3 



+ 



^ + lnf^ 
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Figure 4: Contributions to the scalar stress-energy tensor at order A from 
the conformal counterterm C and the counterterm for the bare cosmological 
constant A. 

-iin=(sj(,,))-^i,i(nM) + ?-^ 

It is an important check of the calculation that ( [L05| ) and ( |106| ) do satisfy 
the D-dimensional conservation law: 

Ppk = -(3 - e)H{ppK + Ppk) ■ (107) 

Although conserved, the potential and kinetic terms still harbor diver- 
gences. We eliminate the nonlocal, "overlapping" divergence by choosing the 
divergent part of the conformal counterterm thus, 

5^ = ^ f— V I ^ + 5i,A . (108) 

We will choose the finite part 5^f^t later so as to achieve maximum simplicity 
in the fully renormalized result. The remaining divergences are proportional 
to so they can be absorbed by a cosmological counterterm. The natural 
choice is, 

i^-^J / Y r^(3-f) c r(3-f) 

SttG ~ \ KhV 8e2 2e2 r(3 - e) 
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r(3-056nt+5A,.n,| , (109) 



where we leave the finite part for later determination. 

Because the conformal and cosmological counterterms cancel all ultravio- 
let divergences we can take e to zero in the renormalized energy density and 
pressure, 



1 



2 



- In' {9) 



7 1, 



2 



1 /r^\ 11 ^ 

IniVt) \ 

^ ^ 9 12 



|-iln(||+6.«,„ 



hi (Si) - 5 



The result can be considerably simplified by making a suitable choice for the 
finite parts of the counterterms, 

7 1 /2u\ 1 vr' 

= + _ n — , 5Af„, = . (112) 

36 12 V/// 18 12 ^ ^ 

The final results are, 

p_ = ^/lln2(l]) + 2^-3_l V^i±l^f]-n-il (113) 

We have verified the conjecture of ref. that the kinetic contributions 
are subdominant to the potential contributions by one factor of ln(f2). This 
suffices to prove that the model violates the Weak Energy Condition on 
cosmological scales. 



p_ + p_ = ^1 _ 1 In (fi) + Iq-^ - 1 V ^fi-"-^ \ . (115) 
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The leading term of this expression is exactly that predicted in equation (49) 
of ref. [Q. The difference is that we have rigorously proved it. We also have 
the sub-dominant corrections and we have a procedure that can be pushed 
to arbitrarily high order. 

Two comments are in order concerning (|TT5|). First, there is a logarithmic 



singularity at In(fi) = Ht = whereas the correct result is that the stress- 
energy tensor vanishes then. The problem seems to be caused by assuming 
ln(f2) ^ e at various places in the reduction. This means we should not 
trust the result for times infinitesimally close to Ht = but it should be 
completely reliable at later times. 

The second comment has to do with the length of time during which we 
can expect the Weak Energy Condition to be violated. The 0^ potential 
obviously has some tendency to force the field back down. One might expect 
this classical effect to eventually balance against the quantum uncertainty 
pressure to result in a constant energy density and pressure.^ We can estimate 
how long the system takes to reach this stage by asking when higher order 
effects become comparable with the two loop result. 

Consider a diagram with 2N external scalar lines. At L loop order the 
number of 0^ interaction vertices is, 

V = L + N-1. (116) 

Each contributes a factor of A, so the diagram goes like A^^^~^. The number 
of internal propagators is, 

P = 2L + N-2. (117) 

Since we are computing a Schwinger diagram, there will be V cancellations 
between + and — variations, which give the ^-function imaginary part of 
the logarithm. However, there are also V temporal integrations, each one of 
which can potentially result in an extra factor of ln(i7). Hence the strongest 
possible effect for the 2A^-point vertex at L loop order is, 

^2"^~A^+^-^(ln(fi))''^"^"\ (118) 



^However, one should be alive to the possibility that stochastic effects result in the 
scalar's further migration up its potential in certain portions of the universe 
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The stress-energy tensor corresponds to = scalar lines so the dominant 
contribution at L loop order is, 

T;:-'""'' - 9,.H' [X\n\n)f'' . (119) 

It follows that perturbation theory breaks down when ln(f2) ~ Since 
A is assumed small we see that the phase of super-acceleration can be made 
to last for an enormous number of e-foldings. Note that all the higher point 
diagrams remain perturbatively weak during this entire period, 

lim ~ . (120) 

It should therefore be valid to use perturbation theory almost up to ln(f2) = 
1/v^. 



Note Added: After the completion of this work we became aware of a 
highly significant paper by Starobinsky and Yokoyama ||2^ in which stochas- 



tic techniques are employed to sum the leading powers of Ht at each order. 
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